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Abstract. The orientation fluctuations of the director of a liquid crystal are 
measured, by a sensitive polarization interferometer, close to the Freedericksz 
transition, which is a second order transition driven by an electric field. Using mean 
field theory, we define the work injected into the system by a change of the electric field 
and we calibrate it using Fluctuation-Dissipation Theorem. We show that the work 
fluctuations satisfy the Transient Fluctuation Theorem. An analytical justification of 
this result is given. The open problems for the out of equilibrium case are finally 
discussed. 

PACS numbers: 05.40.-a,05.70.xz-a 



1. Introduction 

The physics of small systems used in nanotechnology and biology has recently received 
an increased interest. In such systems, fluctuations of the work injected become of the 
order of the mean value, leading to unexpected and undesired effects. For example, 
the instantaneous energy transfer can flow from a cold source to a hot one. The 
probabilities of getting positive and negative injected work are quantitatively related 
in non equilibrium system using Fluctuation Theorem (FT) [I]-jH]- This theorem has 
been both theoretically and experimentally studied in Brownian systems described by 
a Langevin equation [6]-[lI]. The Transient Fluctuation Theorem(TFT) of the injected 
work W T considers the work done on the system in the transient state, i.e. considering 
a time interval of duration r which starts immediately after the external force has been 
applied to the system. For these systems, FT holds for all integration time and all 
fluctuation magnitudes : 



where ks is the Boltzmann constant, T the temperature of the heat bath and p(W T ) is 
the probability density function (PDF) of the injected work W T . S is called symmetry 
function. In a recent paper [15] , some surprising results have been obtained for the TFT 
in a spatially extended system where the fluctuations of the work injected by an electric 
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field into a nematic liquid crystal (LC) have been studied. More precisely, the authors 
find an agreement between experimental results and TFT only for particular values 
of the observation time r. We report in this article experimental results on the same 
system using a different measurement technic. We show that in our experiment, TFT 
holds for all integration time and all fluctuation magnitude. The paper is organized 
as follow. In the section [2, we describe the experimental setup. In the section [3], the 
work and the free energy of the LC system are defined. The experimental results are 
presented in section HJ Finally, we conclude in section [51 

2. Experimental setup 

In our experimental apparatus, we measure the spatially averaged alignment of the LC 
molecules, whose local direction of alignment is defined by the unit vectors n. The 
LC is confined between two parallel glass plates at a distance L. The inner surfaces 
of the confining plates have transparent Indium-Tin-Oxyde (ITO) electrodes, used to 
apply an electric field E. The plates surfaces are coated by a thin layer of polymer 
which is mechanically rubbed into one direction. This treatment of the glass plates 
align the molecules of the LC in a unique direction parallel to the surface (planar 
alignment), i.e. all the molecules near the surfaces have the same director parallel to 
x-axis and n = (1,0,0) (see figured^,). This kind of coating excludes weak anchoring 
effect described in the literature [TTj . The experiments are performed on nematic liquid 
crystal (5CB) produced by Merck in a cell of thickness L = 9 /xm. 

The LC is submitted to an electric field perpendicular to the plates by applying a 
voltage V between the ITO electrodes. In order to avoid electrohydrodynamic effects of 
the motion of the ions invariably present in the liquid crystal, we apply an AC voltage 
at a frequency of f v = 1 kHz (V = y/2V cos(2vr f v t)) [M [19]. When the voltage V 
exceeds a threshold value V c , the planar state becomes unstable and the LC molecules 
try to align parallel to the electric field (fig. Wfr)- This is known as the Freedericksz 

transition which is second order. Above the critical value of the field, the equilibrium 
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structure of the LC depends on the reduced control parameter defined as e = ^% — 1. 
The sample is thermostated at room temperature, i.e. T = 295 K. The stability of the 
temperature is better than 0.5 K over the duration of the experiment. 

The measurement of the alignment of the LC molecules uses the anisotropic optical 
properties of the LC, i.e. the cell is a birefringent plate whose local optical axis is parallel 
to the director n. This optical anisotropy can be precisely estimated by measuring the 
phase shift $ between two linearly polarized beams which cross the cell, one polarized 
along x-axis (ordinary ray) and the other along the y-axis (extraordinary ray). The 
experimental set-up employed is schematically shown in fig. [lb. The beam produced 
by a stabilized He-Ne laser (A = 632.8 nm) passes through the LC cell and it is finally 
back-reflected into the LC cell by a mirror. The lens between the mirror and the LC 
cell is used to compensate the divergency of the beam inside the interferometer. Inside 
the cell, the laser beam is parallel and has a diameter of about 2 mm. The beam 
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Figure 1. a) The geometry of Freedericksz transition: director configuration for 
Vq < V c and director configuration for VqV c . b) Definition of angular displacements 
(0, ip) of the director of one nematic ft. c) Experimental setup. A laser beam passes 
through the LC cell and a polarization interferometer measures the phase shift 4> 
between the ordinary an extraordinary rays |21j . 



is at normal incidence to the cell and linearly polarized at 45° from the x-axis, i.e., 
can be decomposed in an extraordinary beam and in an ordinary one. The optical 
path difference, between the ordinary and extraordinary beams, is measured by a very 
sensitive polarization interferometer [21] . 

The phase shift $ takes into account the angular displacement in xz-plane, 9, and 
in xy-plane, ip (fig. [lb). For Vb larger than the threshold value V c , $ depend only on 9 
and can be written as : 
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with (n , n e ) the two anistotropic refractive indices [H],[l9] an d A = hDq/A is the area 
of the measuring region of diameter Dq in the (x, y) plane. The global variable of our 
interest is ( defined as the spatially averaged alignment of the LC molecules, and more 
precisely : 

2 1 r r rL 
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Close to the critical value V c , using boundary conditions, the space dependance of 9 
along z has the following form : 9 = 9 (x,y,t) sm(irz/L), thus n x = cos9 [HI [191 120] - 
If 6*o remains small, ( takes a simple form in terms of 9q : 
1 
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Figure 2. Typical measurement cycle, a) Dependence of e on time. The perturbation 
Se is switched on at t — 0. b) One realization of the response of £ when the amplitude 
of the control parameter is changed from eo = 0.1985 and 5e = 2 10~ 4 as shown in a). 
The thin line represents the average response (£) done over 2500 experiments. 



After some algebra, we can show that the phase shift is a linear function of £ : 
^ ^ ft n e (n e + n ) \ 

• = M sg— f J (5) 

% = —{n e -n )L (6) 

The phase $, measured by the interferometer, is acquired with a resolution of 24 bits 
at a sampling rate of 16384 Hz and then filtered at 500 Hz in order to suppress the AC 
voltage at fy. The need of such a high acquisition frequency will be explained in the next 
section. The instrumental noise of the apparatus [21] is three orders of magnitude smaller 
than the amplitude <5$ of the fluctuations of $ induced by the thermal fluctuations of 
£. In a recent paper, we have shown that ( is characterized by a mean value (() oc e 
and fluctuations which have a lorentzian spectrum [22l [16] . 

We are interested here in the fluctuations of the work injected into the LC when it 
is driven away from an equilibrium state to another using a small change of the voltage 
(i.e. the control parameter), specifically: 

£ . < < o 

e + 5e t > 

We take care that the system is in equilibrium at t = by applying the value eo for a time 
much larger than the relaxation time (r r ) of the system. A typical measurement cycle 
is depicted in Fig. [2] where we plot e(t) (Fig. EK) and ((t) (Fig. [2b) as functions of time. 
We repeat the same cycle 2500 times to compute, over this ensemble of experiments, 
the probability density function of the work W T injected during a time r. We choose 
the amplitude of 5e sufficiently small such that the mean response of the system to this 
excitation is comparable to the thermal noise amplitude, as can be seen in figure ^jp) 
for e = 0.1985 and 5e = 2 10" 4 . 
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3. Landau equation and definition of the work 



We define now the work injected into this system using the kind of forcing plotted in 
fig. [2^. The free energy is the sum of the electric energy and the potential energy : 
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Elastic free energy is a function of the director n and its spatial derivative [181 [19] and 
electric free energy is f e iect — — l/2eoe a (i?-n) 2 where e a is the dielectric anisotropy of the 
LC (e a is positive for 5CB). As the area of the measuring beam is much larger than the 
correlation length, we can neglect spatial derivative over x and y. Using the sinusoidal 



approximation for 9, 
following form : 



i.e. 



9(x,y,z,t) = 9q(x, y, t) sm(irz/L) [18], free energy has the 
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1 and Kx, K 2 and K3 are the three 
elastic constants of the LC. To simplify the analysis the dependance of #0 in x and y 
can be neglected and the free energy takes the simple form: 

1 , 
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where B 



Air 2 Ki/4L. From eq{TT]we can obtain the equation of momentum 
jAL d9 
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2 dt 

where 7 is a viscosity coefficient and 77 a thermal noise delta-correlated in time, such 
that : 

(ij(t)) = and (r)(ti)r)(t 2 )) = k B T^AL8(t 2 - h) (13) 

where (•) stands for ensemble average. First, we have to calibrate the system, i.e. 
determine the value of the constant B that is a function of the area A which is 
very difficult to precisely determine on the experiment. For this calibration, we use 
Fluctuation-Dissipation Theorem (FDT). We apply a small change of the voltage Vq 
and, so of the control parameter, i.e. e = e + 8e. We separate 9 into the average part 
at eo and a deviation due to 5e : #o(^) = i>o + A(t), where ipQ = 2eo/(ft + e + 1) is the 
stationary solution of eqH2]at e = eo- If the response is linear, the average value of A 
is of the order of 5e. In this limit, the equation of momentum can be rewritten : 

^^ = B[2e -( K + e + ml] A+ 2B5e^Jl-^] 

The external torque, and so the conjugate variable to #o> is equal to 
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We define the integrated linear response function using an heaviside for 5e : 
p/ \ (A(r)) 1 

R{r) = -p = 7 -iV^e (16) 

L ext 2B ip M - ^ 

where _R,5 e is the integrated response function of 6*o to <5e. i? is related by FDT to the 
autocorrelation function, Ce (r) = (59 (t + r) 5#o(£)), of the fluctuations 59 of 6> , that 
is: 

R{r) = ^-{CeM-Ce Q {r)) (17) 
k B T 

The global variable measured by the interferometer and defined in eq. [3] is ( = 
1/Aff 9Q<lxdy ~ + <5#o- Thus the mean value of C is Co = ipo an d the fluctuations 
of C can be related to the fluctuations of 9 : 5( = 2ip 59o. We measure the 
autocorrelation function of S(, = 4^qC*6) . The response function of ( to 5e is also 
related to R$ e '■ R(,5e = ZipoRst- Using fluctuation-dissipation theorem, we can measure 
the constant B : 

-j^Rcse = j^f (Cf (0) - C c (r)) (18) 

where B' = B (1-^/2). 

We define now the work, W T , injected into the system when the control parameter 
e is suddenly changed from e to e = e + Se at t = (eqJTJ). Using equation [T2J W T can 
be written as : 

5e(t')^(29 - 9 3 )dt> (19) 

we consider that the external torque is B5e{t') (29 — 9^). This expression can be 
integrated into : 

W T = B5e ( 9l{r) - 6> 2 (0) - d MzJM ] , (20 ) 



Using the assumption of the limit of small 9q and neglecting its dependance in the x 
and y directions, the work has a simple form in terms of ( = 9q : 

W T = B5e U(t) - C(0) - C2(r) ~ C2(0) ^) (21) 

Notice that if, eo is not too large, the quadratic term in this expression can be neglected, 
but it is not the case in general. The work W T is of course a fluctuating quantity because 
of the fluctuation of £. 



4. Fluctuations of the work injected into the system 

4-1. Experimental results 

The first step consists in calibrating the system. Indeed as we have already mentioned 
the parameter B' appearing in eqlTS] depends on A, i.e. the effective area of the laser 
beam inside the cell that is very difficult to estimate with a good accuracy. In order to 
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Figure 3. Fluctuation-dissipation theorem. g e (r) is represented as a function 
of C<5f(r). Continuous line represents linear fits. The measurement has been done at 
e = 0.1985 with Se = 2 10~ 4 



calibrate the system we use FDT. In this case we use a sufficiently small perturbation 
Se to insure that the response is linear. Thus we obtain the integrated linear response of 
the system : R^s e (r). Then we measure the correlation Q. The function R^s e {r)/B' is 
a linear function of C, as can be seen in figure [3] and the slope seems to be independent 
of (£). Using FDT, the slope is equal to B'/ksT; by this method, we can calibrate 
the constant B' . We measure B' = 4.03 10~ 2 ± 2 10~ 3 . Using in the equation for 
B' = (1 - ^/2)(n 2 K 1 A)/(4:L), the measured values of B' and tp , K x = 6.4 10" 12 N 
(value for the LC 5CB) and L = 9 fim we estimate A = 2.4 10~ 6 m 2 . This value is 
consistent with a, D = 1.75 mm which is very close to the estimated diameter of the 
laser beam inside the cell which is about 2 mm. This justify our approach. 

The system is at equilibrium at t = s and we are now interested in the work 
done by the change of Se during a time r, W T . We need a high acquisition frequency to 
reduce experimental errors in the determination of the first point of each cycle t = 0. 
The probability density function of W T are plotted in figure H] for different values of r. 
We find that the PDFs are Gaussian for all values of r and the average value of W T 
is equal to a few /c#T. We also notice that the probability of having negative values 
decreases when r is increased. The symmetry function is plotted in fig. Hb. It is a linear 
function of W T for any r, that is S(W T ) = H(t)W t . Within experimental errors, we 
measure the slope E(r) = 1. Thus, for our experimental system, the TFT is verified for 
any time r. 

To test the importance of the quartic terms in the expression of the work, we have 
done the same experiments using a higher value of the control parameter eo = 1.039. 
The perturbation remains small Se = 4 10~ 4 (see figure [5]). The quadratic terms in the 
expression of the work cannot be neglected in this situation. We find that the PDFs 
of the work are also Gaussian for all the integration times. This is normal because the 
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Figure 4. Work fluctuations a) PDFs of the work W T injected into the system at 
e-o = 0.1985 by a perturbation 5e = 2 10 -4 , for various r : 10 ms (o), 20 ms (□), 50 ms 
(o) and 100 ms (+). Continuous lines are Gaussian fits, b) Symmetry functions S(W T ) 
computed using the PDF of a) for the same values of r. The straight continuous line 
is a line with slope 1 for all values of t. 



perturbation is such small that we are in the realm of the applicability of the FDT. 
The symmetry function is plotted in fig. [5b. It is a linear function of W T for any r. 
Within experimental errors, we measure the slope equal to 1. Thus, for our experimental 
system, the TFT is verified for any time r. 

4-2. Comparison with theory 

These results can be justified for small eo and Se. Decomposing W T into the sum of a 
mean part (W T ) and a fluctuating one SW T , i.e. W T = (W T ) + 5W T . The distributions 
of the injected work are Gaussian for all values of r, so the symmetry function takes a 
simple form : 

S(W T ) = ^^W T (22) 
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Figure 5. Work fluctuations a) PDFs of the work W T injected into the system at 
e = 1.039 by a perturbation Se = 4 10~ 4 , for various r : 10 ms (o), 20 ms (□), 50 ms 
(+). Continuous lines are Gaussian fits, b) Symmetry functions S(W T ) computed 
using the PDF of a) for the same values of r. The dashed line is a line with slope 1 
for all values of r. 
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where is the variance of the fluctuation of the injected work. The mean value of the 
work is for small €q and small 5e : 



Using FDT, we obtain that 2k B T(W T ) = a^, thus FT is satisfied for all integration 
times and all fluctuation magnitudes as can be seen in fig. lib). 

5. Conclusion 

In conclusion, we have defined in a LC cell close to Freedericksz transition the work 
injected into the system by the change of 5e during a time r using the assumption of 
small angle 8q an d neglecting its dependance in x and y directions. We have shown that 
the TFT holds for the work injected in a spatially extended system for all integration 
times and all fluctuation magnitudes. We have tested this result for another thickness 
of the LC cell and other values of e. It is interesting to notice that this experimental 
test is the first for a spatially extended system in the presence of a non linear potential. 
This is in contrast to the result of ref [15]. The reason of this discrepancy is that they 
neglect in the definition of the work (eq I2Tj) the square term in ( 2 . This cannot be done 
when a large e is used as in the case of ref. [15] . 

Our results are a preliminary study for future investigation on unsolved problems 
such as the applicability of Transient Fluctuation Theorem in an out-of-equilibrium 
system such as an aging system or a rapid quench. 
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